Experimental quantum process tomography of non trace-preserving maps 
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The ability of fully reconstructing quantum maps is a fundamental task of quantum information, 
in particular when coupling with the environment and experimental imperfections of devices are 
taken into account. In this context we carry out a quantum process tomography (QPT) approach 
for a set of non trace-preserving maps. We introduce an operator V to characterize the state 
dependent probability of success for the process under investigation. We also evaluate the result of 
approximating the process with a trace-preserving one. 



P 



in 



Pout 




FIG. 1. (Color online) Scheme of a generic quantum process 
S. 



I. INTRODUCTION 

The complete characterization of quantum devices rep- 
resents one of the fundamental tasks of quantum infor- 
mation science. The characterization of single- and two- 
qubit devices is particularly important since single-qubit 
and two-qubit controlled-NOT gates are the two building 
blocks of a quantum computer Furthermore, identi- 
fying an unknown quantum process acting on a quantum 
system is another key task for quantum dynamics con- 
trol, in particular in presence of decoherence [HQ. In 
this context any quantum process £ can be described by 
a linear map [l| acting on density matrices p associated 
to a Hilbert space T~L which transforms an input state pi n 
into an output state p out (Fig. [I]): 
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Pout — £(Pin)- 



(i) 



The complete characterization of such a process can be 
realized by reconstructing the corresponding map £ . The 
action induced by a black box may be represented by a 
process matrix x which is experimentally reconstructed 
by quantum process tomography (QPT) [3-8]. So far, 
several QPT experiments have been performed for trace- 
preservingprocesses, such as single-qubit transmission 
channels [J [To[ , optimal transpose map \H\ , gates for 
ensembles of two-qubit systems in NMR [12j , a two-qubit 
quantum-state filter [l3| , a universal two-qubit gate Q j 
Control-NOT (CNOT) and Control-Z (CZ) gates for pho- 
tons [IM3- 

Recently theoretical and experimental analyses of non 
trace-preserving processes have been carried out. Kiesel 
et al. evaluated the role of prior knowledge of the intrin- 
sic feature of the experimental setup in order to obtain 



physical and easily understandable parameters for char- 
acterizing the gate and estimating its performance [l7| . 
Furthermore quantum process tomography in presence of 
decoherence has been analyzed for a fast identification of 
the main decoherence mechanisms associated to an ex- 
perimental map [2j. 

Here we address the characterization of non trace- 
preserving maps, focusing on the evaluation of an opera- 
tor representing the success probability of the process. 
In particular we carry out a quantum process tomography 
(QPT) approach for a set of non trace-preserving maps. 
Then, we discuss possible errors occurring in presence of 
inappropriate approximations. 

The paper is organized as follows. In Section II a brief 
review of the main theoretical aspects of QPT and of the 
process fidelity, both for trace-preserving and non trace- 
preserving maps, is presented. In Section III we report 
an example of QPT of a non trace-preserving process, 
corresponding to the transformation induced by a par- 
tially transmitting polarizing beam splitter. The QPT 
experimental realization and results are then presented 
together with a brief discussion on possible wrong ap- 
proaches to the problem, when a non trace-preserving 
process is approximated with a trace-preserving one. Fi- 
nally, the conclusions are given in Section IV. 



II. QUANTUM PROCESS TOMOGRAPHY 

Consider an unknown quantum process, i.e. a black 
box, acting on a physical quantum system described by 
a density matrix p associated to a d-dimensional Hilbert 
space H. A complete characterization of the process may 
be obtained by the Kraus representation of quantum op- 
erations in an open system [3|, 0, EH • A generic map £ 
acting on a generic state p can be expressed by the Kraus 
representation: [l| 



£{p) = Y,E iP E\ 



(2) 
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where Ei are operators acting on the system and satisfy- 
ing the relation^ J2i E\Ei < I. If £ is a trace-preserving 
process the completeness relation J2 i E\Ei = I holds. 

The quantum process tomography of £ consists of the 
experimental reconstruction of the operators {Ei}. In 
order to relate each operator E{ with measurable param- 
eters it is convenient to use & fixed basis of operators {Ai} 
such that: 



Ei — ^ ^ dim-A-m . 



(3) 



By substituting this expression in ([2]), the map reads as 
follows 



(4) 



where Xmn = J2i a>imO>i n - By construction, the matrix xs 
with elements Xmn is hermitian and semidefinite positive. 

To experimentally reconstruct each element Xmn we 
prepare d 2 input states pk forming a basis for the Hilbert 
space ofdxd matrices. The output states can be written 



(5) 



where the coefficients Xkj are experimentally obtained by 
characterizing £{pk) and expressing it in the {pk} basis. 
By defining the coefficients PJ^ 71 such that 

A mPj Ai = J2PTk n Pk, (6) 

k 

it is easy to obtain a relation between Xkj and Xmn 0- 

^2 Pjk n Xmn = X jk . (7) 
m,n 

In order to find the matrix xs which completely describes 
the process £ , we need to operate a matrix inversion 
is this generalized inverse matrix (i.e. 
)m 5 qn ), the elements of xs read: 



of/3$". IfrJT 

l^jk ' jk^jk ~ °1 
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(8) 



For a non trace-preserving map, it is important to con- 
sider not only the transformation acting on a generic in- 
put state, but also the probability of success of the map. 
For a given input state p, the probability of obtaining an 
output state from the black box is given by 



Tv[£(p)}=Tr 



%A m pA^ n 



Tr[Pp] , (9) 



1 EjEi < I means that the eigenvalues of the hermitian oper- 
ator E^Ei — I are not positive. 



where V is a semidefinite positive hermitian operator de- 
fined as: 



V 



mn 



XmnA^Am < I . 



(10) 



Let's write V in its diagonal form, V — J2iPi\Pi)(Pi\, 
where \pi) are the eigenstates and < Pi < 1 the corre- 
sponding eigenvalues. Different cases may occur: 

i) pi = 1 Vi, i.e. V = I for a trace-preserving process; 

ii) pi = p < 1 \/i (V is proportional to the identity 
operator) for a non trace-preserving process with 
state independent success probability; 

hi) there is at least one eigenvalue pi different from the 
others in the case of a non trace-preserving process 
with state dependent success probability. 

The eigenvectors of V coincide with the "probability of 
success" eigenstates of the transformation. 

We now describe how to compare two quantum pro- 
cesses. It is well known that a quantum state can be 
completely determined by a tomographic reconstruction 
[IH and compared with the expected theoretical state 
by a variety of measures, such as quantum state fidelity 
jl9j . Similarly, we know that the process matrix xs gives 
a convenient way of representing a general operation £. A 
closely related but more abstract representation is pro- 
vided by the Jamiolkowski isomorphism [20|, which re- 
lates a quantum operation £ to a quantum state, pg: 



PS = (I®E)\$)($\, 



(11) 



where |3>) = ^= J2j \ 1S a maximally entangled state 
associated to the d- dimensional system with another 
copy of itself, and {\j)} is an orthonormal basis set. If 
£ is a trace-preserving process, then the quantum state 
ps is normalized, Tr [ps] = 1. In this way, by associating 
a quantum process to a quantum state, for two trace- 
preserving processes £ and a Process Fidelity A has 
been defined as follows I2ll-[24l 



A(£,g)=F(p e ,pg) 



(12) 



where T is the quantum state fidelity T = 
Tr \\J yJ~pePQ\fpe\ [121 • ^ * s easv t° demonstrate that, 
by choosing the set A m = |v^|^)0|| as Kraus oper- 
ators, we have ps = Xs, an d? m general, ^(ps^pg) — 
J-(xs, Xg) if an y complete set of operators A! m satisfying 
Tr[A^A^] = ddmn is used (8 mn is the Kronecker delta). 
Thus, if we want to compare an experimental map x with 
the expected one Xid, the process fidelity is 



Tr 



Vx XidVx 



(13) 



The last expression gives the fidelity of density matrices 
with unit trace. However, if x represents a non trace- 
preserving process, i.e. Tr[x] = ^TrfP] < 1, the process 
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fidelity definition is generalized as follows [17j . Let Xid be 
the ideal matrix associated to a non trace-preserving pro- 
cess in the Kraus representation and x the experimental 
one. The fidelity for such a process is written as 



A(x,Xtd) 



Tr [vVx WX 



Tr [ X ] Tr [ Xid ] 



(14) 



Note that the physical meaning of this expression is the 
same of (fT3|) : indeed we can express it as 



A( X , Xjd ) = A( X ',x^)=Tr 



where \' 



(15) 



are well defined phys- 



T^I and X'id ~ T^~ d 

ical states (Tr[x'] = Tr[x- rf ] = 1) which, however, do not 
correspond to any meaningful quantum operation, since 
the probability of success of the corresponding processes 
will be larger than 1 for some input states (i.e. the cor- 
responding operators V will have at least one eigenvalue 
larger than 1). 

It is interesting to highlight that the process fidelity de- 
fined in (p3j) does not distinguish between two processes £ 
and Q if £ = where a is a constant, i.e. two processes 
are indistinguishable if they differ only for a global loss, 
as it often occurs in the experimental implementations of 
photonic quantum systems. 



III. QPT OF A PARTIALLY TRANSMITTING 
POLARIZING BEAM SPLITTER 

Now we analyze a simple example of the quantum pro- 
cess tomography of a non trace-preserving, state depen- 
dent map, acting on a single polarization qubit (d = 2). 
Consider a partially transmitting polarizing beam split- 
ter (PPBS) with trasmittivities Tu and Ty at the hori- 
zontal and vertical polarization, respectively. Following 
the Kraus approach, in which we consider Aj as the Pauli 
operators, we report the analytical expression of the pro- 
cess matrix xppbs- 

In general, if we inject a photon with arbitrary polar- 
ization state into the PPBS the output state will be: 



a\H) + 0\V) -> a^T H \H) + Py/T v \V), 



(16) 



where a, f3 G C and |a| 2 + |/3| 2 = 1. Clearly the probabil- 
ity of success of this transformation is state dependent. 
Let's write the process matrix associated to this map. 
According to QPT calculations, we fix the Pauli matri- 
ces (7^, i = 0, ...,3, in a bidimensional Hilbert space as 
the basis Ai in the Kraus sum (satisfying the normal- 
ization condition Tr[A m At ] = dS mn ), and choose the set 
{pk} of the states to be measured, obtaining the following 
matrix: 



Xppbs 



' 4 




Th-T v 
4 



Th-T v \ 
4 » 





_0 



(17) 



>7 



Obviously, the explicit form of Xppbs does not depend on 
the chosen set but only on the fixed basis Ai in the 

Kraus representation. Let's now write the explicit form 
of the operator V for the PPBS. By using the x matrix 
given in (fTT]h we obtain 



T^ppbs 



T H 
Ty 



(18) 



This operator is proportional to the identity only when 

T H = Ty. 



A. Experimental QPT of a PPBS 

In this subsection we report the experimental realiza- 
tion of QPT for a partially transmitting polarizing beam 
splitter. In the experimental setup shown in Fig. [2] the 
PPBS is implemented by a closed- loop scheme, similar 
to the one used in [25|, [26|, operating with two half- 
waveplates (HWP). A diagonally polarized light beam is 
splitted by a polarizing beam splitter (PBS) in two beams 
with equal intensity and orthogonal polarizations. Pre- 
cisely, the horizontal (H) and vertical (V) components 
travel along two parallel directions inside the interferom- 
eter, counterclockwise and clockwise, respectively. One 
half-waveplate intercepts the H beam, while the other 
intercepts the V beam; by rotating the waveplates it is 
possible to vary the value of Ty with respect to T# . 
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FIG. 2. (Color online) Experimental setup used for the QPT 
of a partially transmitting polarizing beam splitter. The 
PPBS is implemented by a displaced Sagnac interferometer 
and two half- waveplates. The measurements are performed 
with a standard polarization analysis setup. 

The photons injected in this interferometric setup are 
generated by a spontaneous paramentric down conversion 
source realized with a nonlinear crystal cut for type II non 
collinear phase matching [27j. The crystal is pumped 
by a CW diode laser and pairs of degenerate photons 
are produced at wavelength A = 806nra. One of the 
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photon is used as a trigger, while the other is delivered 
to the PPBS setup. We prepared six different input 



1.025 



function 
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a, 6=1 



nab- ^2 (V^ml^a)' 
m,n=0 

{4>a\d-„\lpb)Xmn(i) 



(19) 



where n a b are the measured coincidence counts for the 
ath input and the 6th output, \<p a ) and indicates the 
input and the output state respectively, and a m are the 
Pauli operators. Since we are not interested into overall 
losses affecting the transformation (even the adopted fi- 
delity is independent of global losses) we normalize the 
experimental Xexp matrix such that the maximum eigen- 
value of V is 1. We determined the fidelity between the 



FIG. 3. (Color online) Measurement of fidelity as a function 
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FIG. 4. (Color online) Real part of ideal and experimental 
process matrices, Xid and \exp-, for (a-b) T = 0.879 and (c-d) 
r = 0.255. The imaginary part are negligible. 

states, | if), \V), \D), \A), \R), \L) associated to horizon- 
tal, vertical, diagonal, anti-diagonal, right-handed and 
left-handed polarization respectively, and measured the 
six output components for each input with a standard 
polarization analysis setup. We repeated this procedure 
for different values of the ratio T = Ty/Tn and, for each 
value of T, we reconstructed the experimental \ matrix 
of the process. We then performed an optimization of 
the process matrix following a maximum likelihood ap- 
proach [l5|, [28| ; in particular we minimized the following 




FIG. 5. (Color online) Probability operator eigenvalues as 
a function of the ratio Y — Tv/Th- Solid lines represent 
expected behaviour. 

experimental map and the ideal one for several values 
of T, as shown in Fig. 03 We observe that the pro- 
cess fidelity approaches unity for each value of T, and in 
general, we have F > 96% with a good agreement be- 
tween the experimental data and the theory. In Fig. [H 
two examples of ideal and experimental process matrices, 
corresponding to T = 0.879 and T = 0.255, are shown. 

We also estimated the probability operator V: the be- 
haviour of its eigenvalues Ai and A2 as a function of T is 
shown in Fig. [5l We observe that Ai = 1 for each value 
of T (by construction), while the other eigenvalue, A2, 
shows a decreasing behaviour as the ratio between the 
trasmittivities decreases, as expected from (fl8|) . Again, 
a very good agreement between experimental data and 
theory is obtained. 

B. Trace-preserving approximation 

The method above described can be usefully adopted 
even when the process under investigation is ideally 
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FIG. 6. (Color online) Experimental fidelities calculated im- 
posing the constraint n Xmncr n crm = I (red open circles). 
Fidelities obtained with the previous method are also reported 
(black filled triangles). 
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(fT7|) for each value of T. The results are shown in Fig. [U 
As expected, this method gives results similar to those 
obtained in Section IIIA for T ~ 1, while the fidelities 
values are different as T decreases. In particular, the fi- 
delities calculated by imposing the constraint decrease as 
T goes to zero. It is evident that constraining the pro- 
cess to be trace-preserving does not allow to correctly 
reconstruct the associated map. 

A further scenario where probability of success must 
be taken into account may arise when measurements are 
performed in post-selection. The reconstruction of the 
output state density matrices (which obviously are nor- 
malized physical states) for several input state, leads to 
a trace-preserving process. Even in this case we evalu- 
ated the fidelities between the resulting process matrix 
and the ideal one obtaining the results shown in Fig. 
[7l As in the previous case the fidelity decreases as T 
goes to zero. Note that this approach is not correct 
even from a theoretical point of view: the process ma- 
trix xs obtained by normalizing the output states could 
be non-physical (i.e. it could have negative eigenvalues) 
and its expression depends on the chosen set of input 
states. This is due to the fact that normalization implies 
the process to be no longer a linear map and equation (J4j) 
is not valid anymore. In general, the output state nor- 
malization produces wrong process matrices for any non 
trace-preserving operation with state dependent success 
probability. 



IV. CONCLUSIONS 



FIG. 7. (Color online) Experimental fidelities calculated us- 
ing the post-selective approach (red open circles). Fidelities 
obtained with the correct method are also reported (black 
filled triangles). 



trace-preserving. In fact, when a quantum process 
tomography is practically implemented, any interac- 
tion with the environment as well as experimental im- 
perfections may cause the process to be non trace- 
preserving. In practice, to approximate the process 
as a trace-preserving one corresponds to minimize the 
likelihood function ([T9]) with the additional constraint 
V J2m,n Xmncrncrm I. In this way we are imposing 
the probability of success to be independent of the input 
state. We carried out the fit) minimization by taking 
into account the constraint^ and evaluated the process 
fidelity between the obtained Xexp and the ideal matrix 



2 We used the function Minimize [{f , cons}, t] of the 
MATHEMATICA©5 program that allows to numerically mini- 



mizes f(t) subject to the constraints cons. V - 
constraint imposes the normalization Tr[xeccp] 



Note that the 



A review on quantum process tomography of non trace- 
preserving maps has been reported. The experimental 
implementation of a simple non trace-preserving, state 
dependent process, i.e. the transformation induced by a 
partially polarizing beam splitter, provided process fideli- 
ties larger than 96% for any value of the ratio between 
the transmittivities T. Particular attention has been ad- 
dressed to the state dependence property of the process 
through evaluation of the operator V ([T0|) . This oper- 
ator has been calculated and measured in the case of a 
PPBS and its eigenvalues resulted to be different from 
unity (see ([18]) ). as expected for a non trace-preserving 
process. In order to stress the validity of the method 
a brief discussion about possible wrong approaches has 
been presented together with the explicit calculation of 
the PPBS process fidelities. The obtained results clearly 
show that approximation of a non trace-preserving, state 
dependent process with a trace-preserving one does not 
allow a correct reconstruction of the real process map. 

QPT of non trace preserving processes are relevant for 
linear optical logic gates with success probability < 1. 
Indeed, typically it is just assumed that the success prob- 
ability of such gates is uniform across input states and 
hence it is crucial to check the validity of this assumption 
for any application. For example, it can be interesting to 
investigate whether losses in the planar integrated waveg- 



uide chips currently being used [29| could affect different 
input states differently. 
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